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ABSTRACT 

Gamma ray bursts (GRBs) have recently attracted much attention as a possible way 
to extend the Hubble diagram to very high redshift. To this aim, the luminosity (or 
isotropic emitted energy) of a GRB at redshift z must be evaluated from a correlation 
Q ■ with a distance independent quantity so that one can then solve for the luminosity dis- 

^ ' tance Dl{z) and hence the distance modulus /Lt(z). Averaging over five different two 

C/3 , parameters correlations and using a fiducial cosmological model to calibrate them, 

Schaefer (2007) has compiled a sample of 69 GRBs with measured n(z) which has 
since then been widely used to constrain cosmological parameters. We update here 
. that sample by many aspects. First, we add a recently found correlation for the X- 

^ ' ray afterglow and use a Bayesian inspired fitting method to calibrate the different 

I GRBs correlations known insofar assuming a fiducial ACDM model in agreement with 

0^ ■ the recent WMAP5 data. Averaging over six correlations, we end with a new GRBs 

Hubble diagram comprising 83 objects. We also extensively explore the impact of 
varying the fiducial cosmological model considering how the estimated ^{z) change 
as a function of the {il.M,wo,Wa) parameters of the Chevalher- Polarski- Linder phe- 
' nomenological dark energy equation of state. In order to avoid the need of assuming 

0^ . an a priori cosmological model, we present a new calibration procedure based on a 

■ model independent local regression estimate of fi{z) using the Union SNela sample to 

calibrate the GRBs correlations. This finally gives us a GRBs Hubble diagram made 
- out of 69 GRBs whose estimated distance modulus fi{z) is almost independent on the 

' underlying cosmological model. 

H ; 

Cq ^ Key words: Gamma Rays : bursts - Cosmology : distance scale - Cosmology : cos- 

mological parameters 



1 INTRODUCTION 

The Hubble diagram of Type la Supernovae (SNela) pro- 
vided the first evidence of the cosmic speed up (Perlmutter 
et al. 1997; Riess et al. 1998) opening the way to a flood 
of theoretical and observational efforts devoted to investi- 
gate the nature and the nurture of such a phenomenon. The 
astrophysical data coming from later SNela samples (Riess 
et al. 2007; Astier et al. 2006; Wood-Vasey et al. 2007; ?; 
Davis et al. 2007; Kowalski et al. 2008), Baryonic Acoustic 
Oscillations (Eisenstein et al. 2005; Percival et al. 2007), the 
anisotropy spectrum of cosmic microwave background radi- 
ation (de Bernardis et al. 2000; Bennett et al. 2003; Bennett 
et al. 2003; Dunkley et al. 2008) and the large scale struc- 
ture data from large galaxy redshift surveys (Dodelson et 



al. 2002; Percival et al. 2002; Szalay et al. 2003; Hawkins et 
al. 2003) have strenghtened the then surprising result. As a 
consequence, there is nowadays little doubt that the universe 
is spatially flat, has a low matter content {Qm — 0.3) and 
is undergoing a phase of accelerated expansion. To close the 
energy budget and drive the cosmic speed up, a new fluid 
with negative pressure dubbed dark energy has been invoked 
as the only viable solution. Understanding what dark energy 
really is stands out as one of the most fascinating problems 
of modern cosmology. Roughly speaking, we can divide the 
different proposals in two broad categories. In the flrst class, 
one can include all models considering dark energy as a true 
fluid with unusual properties with proposals ranging from 
the classical cosmological constant A (Carroll et al. 1992; 
Sahni & Starobinski 2000) to scalar fields running down their 
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self interaction potential (Peebles & Rathra 2003; Padman- 
abhan 2003; Copcland et al. 2006). On the other hand, one 
may give off General Relativity resorting to modified grav- 
ity theories such as f{R) models (see, e.g., Capozziello et al. 
2008, Sotiriou & Faraoni 2008 and references therein) and 
braneworld cosmologies (Dvali et al. 2000; Lue et al. 2004) 
manifesting themselves as a fictitious dark energy - like fluid. 

Notwithstanding their radically different inspiring phy- 
losophies, all the viable theories share the same ability to 
be in satisfactory agreement with the available data. Such a 
confusing situation is a clear evidence that present dataset 
are unable to discriminate among different models. Such a 
problem comes out as a consequence of the still low redshift 
range probed. The SNela Hubble diagram is an illuminating 
example. As well known, the distance modulus fi{z) scales 
linearly with the logarithm of the luminosity distance dz,(«) 
which depends on the dark energy equation of state through 
a double integration. Discriminating among different models 
asks for extending the Hubble diagram to higher redshifts 
since the expressions for the H{z) and hence fj,{z) are more 
and more different as one goes to higher and higher z values. 

Unfortunately, such a project can not be undertaken 
relying on SNela since, even with excellent space based 
projects such as SNAP (Aldering et al. 2004), it is almost 
impossible going to 2: > 2. Indeed, up to now, most of the 
SNela observed have redshift z < 1.2 with only one SN at 
z = 1.7. On the contrary, Gamma Ray Bursts (GRBs) are 
ideal candidates from this point of view. Indeed, thanks to 
their enormous energy release, they have been observed up 
to redshift z > 5, the most distant one being GRB 080913 
with z = 6.7 (Greiner et al. 2008). as high as z ~ 6.5. There 
is, however, a fierce problem to overcome. Actually, GRBs 
arc not standard candles so that it is not immediately pos- 
sible to build their Hubble diagram. To this aim, one has to 
look for a correlation between a distance dependent quantity 
and a directly observable property. Many of these correla- 
tions have been suggested in the last few years so that the 
route towards making GRBs standardizeable candles (in the 
same way as SNela) is now open. 

Actually, various attempts to build a GRBs Hubble di- 
agram have yet been made. Bradley (2003) first derived the 
luminosity distances of 9 GRBs with known redshifts by us- 
ing two GRBs luminosity relations thus providing the first 
GRBs Hubble diagram. Dai et al. (2004) and Xu et al. (2005) 
have then proposed a methodology to constrain cosmological 
parameters using GRBs and apply it to preliminary samples 
(made out of 12 and 17 objects respectively) relying on the 
Ghirlanda relation. A Bayesian based approach was instead 
proposed by Firmani et al. (2005), later improved by the 
same authors (Firmani at al. 2006) advocating the joint use 
of GRBs and SNela. 

As a Art summary of these efforts, Schaefer (2007) has 
compiled a catalog of 69 GRBs with measured properties 
entering the five two parameters correlation then available. 
Estimating n{z) from each single correlation and performing 
a weighted average, Schaefer have finally built the first GRBs 
Hubble diagram which has later been used in different cos- 
mological analyses (Wright 2007; Wang et al. 2007; Li et al. 
2008; Daly et al. 2008; Capozziello & Izzo 2008). We update 
here the Schaefer work by many aspects. First, we add the 
recently found Lx Ta correlation (Dainotti et al. 2008) with 
Lx the X - ray luminosity at the time Ta and Ta a timescale 



characterizing the late afterglow decay. The use of this new 
correlation allows us to both increase the GRBs sample and 
reduce the uncertainty on fj,{z). We then recalibrate all the 
six correlations considered using a fiducial ACDM cosmo- 
logical model in agreement with the WMAP5 (Dunkloy et 
al. 2008) data and a Bayesian motivated fitting technique 
(D' Agostini 2005). We are now able to build a new GRBs 
Hubble diagram whose dependence on the background cos- 
mological model used for the calibration is explored. Finally, 
we present a novel method to escape the circularity prob- 
lem, i.e. to calibrate the GRBs empirical laws without as- 
suming any a priori cosmological model. This gives us then 
the possibility to build up a new GRBs Hubble diagram less 
affected by systematic uncertainties related to the unknown 
cosmological scenario. 

The plan of the paper is as follows. In Sect. 2, we briefly 
review the different GRBs 2D correlations known insofar 
and present the Bayesian motivated method wo will use in 
the following to calibrate them. Such a calibration is per- 
formed in Sect. 3.1 using a fiducial ACDM model with pa- 
rameters set according to the recent WMAP5 results. These 
calibrated relations are then used in Sect. 3.2 to build a 
Hubble diagram comprising 83 objects, while Sect. 3.3 ex- 
plores what is the impact of varying the cosmological model 
used in the calibration on the estimated distance modulus. 
Sect. 4.1 presents the local regression method showing how 
it is possible to use the Union SNela sample to recover the 
value of ii{z) in a cosmological model independent way. This 
method is then used to calibrate the 2D correlations for the 
GRBs with z < 1.4 in Sect. 4.2 so that a newly calibrated 
GRBs Hubble diagram may be built in Sect. 4.3 where we 
also compare this latter to the fiducial one. The impact of 
different kind of systcmatics on the derived Hubble diagrams 
is investigated in Sect. 5. Having thus checked the reliability 
of the GRBs HDs thus obtained, we present in Sect. 6 a first 
cosmological analysis and discusses the importance of GRBs 
as a cosmological tool presenting a forecast of the precision 
affordable on the matter density parameter, Qm, and the 
dark energy equation of state. Finally, Sect. 7 is devoted to 
the conclusions, while in the Appendix we tabulate the es- 
timated distance modulus from the two different methods 
and make some few comments on possible outliers. 

2 TWO PARAMETERS CORRELATIONS 

In order to be as general as possible, let us denote with x a 
distance independent property and let y = kd\{z) a given 
quantity with k a redshift independent constant and dz,(«) 
the luminosity distance. Let us suppose that we have been 
able to find a correlation between them such as : 

log y — k + 2 log dL{z) = a -\~ 6 log x . 

Should the relation be calibrated, i.e. the parameters (a, b) 
be known, we can then use it to compute the distance mod- 
ulus /^(z) of an object at redshift z simply as : 

ii{z) = 25 + 5log dL{z) 

= 25 + {5/2){\ogy - k) 

= 25 + {5/2){a + blogx-k) 
where all the quantities are now known. 
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Notwithstanding its conceptual simplicity, performing 
such a program for GRBs is far to be trivial. First, one has 
to find such a two parameters (hereafter, 2D) correlation, i.e. 
decide what are the quantities {x,y). As we will see, there 
arc different choices which seems to work satisfactorily well, 
but the intrinsic scatter around the best fit line is usually 
not negligible. Moreover, it is still not clear what is the un- 
derlying physical mechanism originating the empirical 2D 
correlations found. As a second more important problem, 
calibrating such relations is a still open problem. Indeed, 
since local GRBs, i.e. GRBs with z < 0.01, are not avail- 
able, the only one being GRB 980425 (Galama et al. 1998), 
one has typically assume an a priori cosmological model to 
compute dhiz) and hence y so that the calibration turns 
out to be model dependent thus leading to the well known 
circularity problem when one aims at using GRBs as cosmo- 
logical probes. As a third problem, the calibration procedure 
is performed using a variety of statistical methods so that 
the results from the use of different correlations can not be 
straightforwardly compared. 

In the following, we first briefly review the known 2D 
correlations and then describe a Bayesian motivated fitting 
procedure to determine the calibration parameters (a, h) and 
the intrinsic scatter aint once a GRB sample is available. 



2.1 2D empirical laws 

The 2D correlations we will discuss here typically relate a 
GRB observable with the isotropic absolute luminosity L or 
the coUimation corrected energy E^, with a single case cor- 
relating with the X - ray afterglow luminosity. On the other 
hand, the observable properties cover a wide range of GRBs 
features including both quantities related to the energy spec- 
trum and the light curve. In the first class, we find the peak 
energy Ep, which is the energy at whith the uFi, spectrum 
is brightest, and the time lag Tiag, which measures the time 
offset between the arrival of the low and high energy pho- 
tons. To the light curve shape are instead related the rise 
time trt, which is the shortest time over which the light 
curve rises by half the peak flux of the burst, and the vari- 
ability V , which quantifies the smoothness of the light curve 
itself. Somewhat different is the Lx — Ta correlation, with 
Lx the luminosity at the time Ta and Ta a time scale, both 
referring to the X - ray afterglow decay curve. 

As yet said, neither L nor E^, are directly measur- 
able quantities since they depend on the luminosity distance 
dhiz)- Indeed, it is: 



L = A-Kd\{z)Pboio , 

E^ = 4:nd%{z)SboloFbeam{i- + z) 



(1) 

(2) 



Here, Pboio and the bolometric peak flux and flu- 

encc, respectively, while Fbeain = 1 — cos 6 jet is the beaming 
factor with Ojet the rest frame time of achromatic break 
in the afterglow light curve. Note that Pboio and Sboio are 
computed from the observed GRB energy spectrum <^{E) as 
follows (Schaefer 2007) : 



Pboio = P X 



Si 



10*/(l-|-z) 

/(1+^) 



E<i>{E)dE 



iE:::mE)dE 



(3) 



nolo 



= S X 



ur/(l + z) 



E<!>{E)dE 



(4) 



with P and S the observed peak energy and fluence in units 
of erg/cm^ /s and erg/cm^, respectively, and {Emin,Emax) 
the detection thresholds of the observing instrument. Note 
that the energy spectrum is modelled using a smoothly bro- 
ken power -law (Band et al. 1993) which reads: 

r ^^ae-(2+°)^/'Bp E <[{a- f3)/{2 + a)]Ep 
$(£) = <^ .(5) 

[ BE'^ otherwise 

While L and E^ are bolometric quantities, we will also con- 
sider a recently found 2D correlations involving only X - ray 
observables. In particular, we will therefore be concerned 
with the luminosity in the X-ray. As a general rule, if we 

denote by Lx{t) and Exit) the luminosity and the flux at 
the time t in the afterglow light curve, one can write : 



Lx(t) = 47rdi(z)Fx(t) 



(6) 



where the flux must be corrected (Bloom et al. 2001) as : 

rs„„x/(i+z) 



Fx{t) = m X 



r 

J E 



E^x{E)dE 



E(S>x{E)dE 



(7) 



with (Ernin, Emax) ~ (0.3,10) keV sct by the instrument 
bandpass and ^x{E) the energy spectrum in the X-ray 
band, well approximated as a single power law. Finally, f{t) 
describes the time variation of the flux along the afterglow 
curve and is parametrized as given in Willingale et al. (2007). 
Note that Eq.(6) is the same as Eq.(l) the only difference be- 
ing the integration limits at the numerator. Actually, while 
L is the bolometric luminosity, Lx refers to the X-ray one 
so that we integrate only over this energy range. 

As an important remark, it is worth stressing that all 
the 2D correlations we will use refer to quantities evaluated 
in the GRB rest frame, while the observed spectra and light 
curves are in the Earth orbiting satellite frames. As such, 
they axe affected by time dilation and redshifting which have 
to be corrected for before the calibration procedure. All of 
the 2D correlations available insofar are power - laws so that 
it is useful to work in logarithmic units. As a general rule, 
we will therefore consider a linear rule as : 



log R = a log Q + b 



(8) 



where J? is a distance dependent quantity, while Q is not. 

Different choices for (_R, Q) determine the different 2D em- 
pirical laws. Setting R = E~f and Q = Ep{l + 2:)/300 keV 
(where the factor (1 -I- 2) corrects for redshifting) gives the 
Ghirlanda relation (Ghirlanda et al. 2004; Ghirlanda et al. 
2006). It is worth noting that this is the only relation in- 
volving an energy quantity, the other ones always involving 
a luminosity as R quantity. Indeed, setting R = L, we can 
find 2D correlations for Q — Ep{l + 2;)/300 keV (Schaefer 
2003; Yonekotu et al. 2004), Q = Tiag{l + /O.I s (Nor- 
ris et al. 2000), Q = TBrii- + zy /O.l s (Schaefer 2007) 
and Q ^ V{1 + z)/0.02 (Fenimore & Ramirez- Ruiz 2000; 
Riechart et al. 2001; Schaefer 2007). Note that in aU these 
correlations the factor (1 -|- z) or (1 -|- z)~^ corrects for time 
dilation or redshifting, while the normalization is chosen in 
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order to minimize the correlation between the errors. Fi- 
nally, setting R = Lx{Ta) and Q = Ta/{l + z) gives the 2D 
correlation empirically found for the first time in Dainotti 
et al. (2008) and later confirmed by Ghisellini et al. (2008) 
on a scrnithcorctical basis and by Yamazaki (2008) in the 
framework of his model proposed to explain the observed 
plateau in the X-ray GRBs afterglows. 



2.2 Bayesian fitting procedure 

Eq.(8) is a linear relation which can be fitted to a given 
dataset {Qi,Ri) in order to determine the two calibration 
parameters (o, 6). Actually, the situation is not so simple as 
one may naively expect. Indeed, both the {Q, R) variables 
are affected by measurement uncertainties {<JQ,aR) which 
can not be neglected. Moreover, ctq/Q ~ (tr/R so that the 
simple rule to choose as independent variable in the fit the 
one with the smallest relative error can not be applied here. 
Finally, the 2D correlations we are interested in are still not 
definitively explained by an underlying theoretical model de- 
termining the detailed features of the GRBs explosion and 
afterglow phenomenology. Both on theoretical and observa- 
tional grounds, we do expect a certain amount of intrinsic 
scatter around the best fit line which has to be taken into 
account and determine together with (a, b) by the fitting 
procedure. Different statistical recipes are available to cope 
with these problems with each author having its own pref- 
erences. How the fitting technique employed affects the final 
estimate of the distance modulus for a given GRB from the 
different 2D correlations is not clear so that it is highly de- 
sirable to recalibrate all of them with the same method. 

To this aim, in the following we will resort to a Bayesian 
motivated technique (D' Agostini 2005) thus maximizing the 
likelihood function £(a, 6, aint) = exp [—L{a, b, aint)] with : 

L(a, b, aint) = ^'^^^ i'^int + crii + a^o-Q . ) 

where the sum is over the A/" objects in the sample. Note 

that, actually, this maximization is performed in the two 
parameter space {a, aint) since b may be estimated analyti- 
cally as : 



b = 



Ri — aQi 



2^2 



E: 



+ a-'a, 



2^2 



(10) 



so that we will not consider it anymore as a fit parameter. 

It is worth noting that the Bayesian approach allows to 
find out what is the most likely set of parameters within a 
given theory, but does not tell us whether this model fits 
well or not the data. An easy way to quantitatively estimate 
the goodness of the fit is obtained considering the median 
and root mean square of the best fit residuals, defined as 
5 = Robs — Rfit which we will also compute for the different 
2D correlations we will consider. 

The Bayesian approach used here also allows us to quan- 
tify the uncertainties on the fit parameters. To this aim, for 
a given parameter pi, we first compute the marginalized like- 
lihood jCi{pi) by integrating over the other parameter. The 
median value for the parameter pi is then found by solving : 



/ £i(pi)dpi =2 Ci{pi)dpi . (11) 

The 68% (95%) confidence range {pi,i,pi,h) are then found 

by solving : 

/ jCi{pi)dpi = / Ci{j>i)dpi , (12) 

''Pi,l •'Pi,min 
rPi,h 



Ci{pi)dpi = 



Ci{pi)dpi 



with e = 0.68 (0.95) for the 68% (95%) range respectively 

3 UPDATING THE GRB HUBBLE DIAGRAM 

Following the classical approach, we calibrate the six 2D cor- 
relations described above by first estimating the distance de- 
pendent quantity (namely, E~,,Lov Lx) in a given cosmolog- 
ical model. In particular, to be in agreement with the most 
recent results, we choose a fiducial spatially flat concordance 
ACDM model with = 1 - = 0.291 and h = 0.697 as 
suggested by the analysis of the WMAP5 data (Dunkley et 
al. 2008). The sample of GRBs used is the one compiled by 
Schaefer (2007) for the first five 2D correlations, while the 
Lx - Ta law is calibrated using a subsample of the Willingale 
et al. (2007) catalog made out of 28 GRBs with measured 
redhisft, log Lx{Ta) > 45 and 1 < log [Ta/(1 + «)] < 5 (see 
Dainotti et al. 2008 for the motivation of these cuts). 

3.1 The calibration parameters 

The Bayesian fitting procedure described above then gives us 
the results summarized in Table 1. It is worth stressing that, 
since the calibration parameters have been obtained using 
the same statistical analysis for all the 2D correlations, it is 
now possible to compare the relations avoiding any possible 
bias introduced by the different fitting procedure. 

As a first issue, one can qualitatively judge what is 
the best correlation, i.e. what is the correlation giving the 
smaller residuals or having the smaller intrinsic scatter. 
From an observational point of view, there is not a 2D law 
working significantly better than the other ones since both 
5med and Srms take almost the same values over the full set, 
with only a modest preference for the E^-Ep, L-Tiag and 
Lx - Ta correlations. On a theoretical side, one can argue 
that the 2D law having the smallest intrinsic scatter is the 
better motivated one since the (unknown) underlying phys- 
ical mechanism works in the same way for all the GRBs 
involved. Taken at face values, the maximum likelihood es- 
timates for aint should argue in favour of the - Ep and 
Lx - Ta correlations, but this conclusion becomes meaning- 
less when one takes into account the 68% and 95% confi- 
dence ranges which overlap quite well. 

Finally, we caution the reader that we have implicitely 
assumed that the calibration parameter do not change with 
the rcdshift notwithstanding the fact that z spans a quite 
large range (from z = 0.125 up to z = 6.6). The limited 
rmmber of GRBs makes it not possible to explore in de- 
tail the validity of this usually adopted working hypothesis, 
but we can get a hint by considering whether the residuals 
correlate with the reshift. The values reported in Table 1 



An updated GRBs Hubble diagram 5 



Table 1. Calibration parameters (a, 6) and intrinsic scatter ai„t for the 2D correlations logR = alogQ + b in a fiducial ACDM model. 
Columns are as follows : 1. id of the correlation with the first letter referring to the R variable and the second to the Q quantity; 2. 
number of GRBs used; 3. maximum likelihood parameters; 4, 5, 6. median, root mean square and Spearman rank correlation parameter 
with the redshift of the best fit residuals; 7, 8. median value and 68 and 95% confidence ranges for the parameters (a, ai„t). 

Id A/" {a,b,r7int)ML 5med irms C{5,z) atl^ -2ct {<^int)t.\'^ -2^ 

E-y-Ep 27 (1.38,50.56,0.25) 0.01 0.38 -0.17 -^Tto^l 0.301°; Jg l° ig 



L-Ep 64 (1.24,52.16,0.45) -0.05 0.51 0.07 



4+0.18 +0.36 n y|o+0.07 +0.17 



L-r,ag 38 (-0.80,52.28,0.37) 0.02 0.40 0.40 -O-SOto^t to^^, — -o.07 -0.12 

L-TRT 62 (-0.89,52.48,0.44) 0.04 0.47 0.27 -0.89lg;i| 0.mt.o °^ to,\l 

L-V 51 (1.03,52.49,0.48) -0.09 0.50 0.25 l-O^l^Jg to.5? O-^llao? ^0.13 

Lx-Ta 28 (-0.58,48.09,0.33) -0.12 0.43 -0.21 -0.58tQ;Jg 1^ 0.391°;" IqIo 



shows that only a weak and not significant correlation ex- 
ists thus arguing in favour of the starting hypothesis. This 
is in agreement with what has recently been found by Basi- 
laJios & Perivolaropoulos (2008) which have calibrated the 
first five 2D correlations in Table 1 dividing the GRBs in 
redshift bins. Although they used a different fitting tech- 
nique so that our results can not be directly compared to 
their ones, we nonetheless agree with their conclusion that 
no statistical evidence of a dependence of the (o, h, aint) pa- 
rameters on the redshift is present. 

3.2 Making up the GRBs Hubble diagram 

Once the six 2D correlations have been calibrated, we can 
now use them to compute the GRBs Hubble diagram. In or- 
der to reduce the error and, in a sense, marginalize over the 
possible systematic biases present in each one of the correla- 
tions, we follow Schaefer (2007) and take a weighted average 
of the distance modulus provided by each one of the six 2D 
laws in Table 1. As a preliminary step, let us remember that 
the luminosity distance of a GRB with redshift z may be 
computed as : 



dL{z) = < 



E^{1 + Z) /4n FbeamSbolo 
L/A-KPbolo 

Lxil + zf^+^/A'wFxiTa) 



(14) 



depending on whether a correlation involving E^, L or Lx is 
used*. The uncertainty on dhiz) is then estimated through 

the naive propagation of the measurement errors on the in- 
volved quantities. In particular, the error on the distance 
dependent quantity R is estimated as : 



CT(logii) = VaV2(l0g Q)+fT|„ 



(15) 



* Note that the factor (1 + z)^'^'^^ in the formula involving Lx 
come from the k- correction of the flux Fx{Ta) with Pa being the 
slope of the power - law energy spectrum in the X - ray band. 



and is then added in quadrature to the other terms entering 
Eq.(14) to get the total uncertainty. The distance modulus 
fj,{z) is easily obtained from its definition : 



n{z) = 25 + 5\ogdL{z) 



(16) 



with its uncertainty obtained again by error propagation. 
Following Schaefer (2007), we finally estimate the distance 
modulus for the i - th GRB in the sample at redshift Zi as : 



IJ.{Zi) = 



^Hj{zi)/al. 



with the uncertainty given by : 
-I -1 



(17) 



(18) 



where the sum runs over the 2D empirical laws which can 
be used for the GRB considered. Joining the Willingale et 
al. (2007) and Schaefer (2007) samples and considering that 
17 objects are in common, we end up with a catalogo of 83 
GRBs which we use to build the Hubble diagram plotted 
in Fig. 1 and tabulated in the Appendix^. We will refer in 
the following to this dataset as the fiducial GRBs Hubble 
diagram (hereafter, HD) since it relies on a calibration made 
using a fiducial ACDM model to compute the distances. 

It is worth comparing our fiducial HD to the one de- 
rived by Schaefer (2007) which we will refer to as the Schae- 
fer HD. Actually, we have improved over the Schaefer HD 
by three aspects, namely updating the ACDM model pa- 
rameters, using a Bayesian motivated fitting procedure and 
adding a further 2D correlation. To investigate the impact 
of these differences, we first consider the 69 Schaefer GRBs 
and compute the distance modulus with our new calibration. 



t ASCII tables with the Hubble diagram and all the other 
datasets necessary to repeat the full analysis presented in this 
paper together with the Mathematica codes used are available on 
request to V.F. Cardone (winnyenodracSgmail.com). 
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but not including the values coming from the Lx - Ta corre- 
lation. We then plot the in left panel of Fig. 2 the old vs the 
new estimates of the distance modulus. As can be seen, there 
is a clear one-to-one correspondence with slope perfectly 
compatible with 1 and no systematic offset. This can also 
better appreciated looking at the right panel showing that 
l-i'nem / Hold IS closB to 1 withiu the 5% with no correlation 
with the redshift. To be more quantitative, averaging over 
the 69 GRBs, we find fi„ew/f^oid ~ l.OOiO.Ol with an root 
mean square value {i-Lnem/fJ'oid)rm3 = 1.002. We may there- 
fore safely conclude that the new calibration procedure have 
not affected the results. 

A cautionary remark on the errors is in order here. As 
explained above, the uncertainty on fi comes from the prop- 
agation of the errors starting from the ones on the distance 



dependent quantities E.y, L or Lx- Comparing Eq.(15) with 
the similar formulae in Schaefer (2007), one can see that we 
have not included in the propagation the errors on the cali- 
bration coefficient {a,b). This is a direct consequence of the 
different statistical phylosophy underlying the fitting proce- 
dure employed (Bayesian maximum likelihood here vs the 
frequentist minimization in Schaefer). We have however 
checked whether this bias somewhat the error estimates. Av- 
eraging over the 69 GRBs, we find anew/iyoid = 1.01±0.06 
with a root mean square value {<J„ew /o'oid)rms = 1.01 and 
no correlation with the redshift. These numbers make us 
then confident that no under or overestimation of the errors 
is introduced by our different approach. 

Having added the new Lx - Ta correlation makes it pos- 
sible to infer the distance modulus for further 14 objects thus 
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Figure 3. Root mean square of the percentage deviation of the distance modulus /i(2) (left) and the luminosity distance d^iz) from the 
fiducial ACDM values caused by changing the {wo,Wa) parameters for Qm = 0.30. Similar results are obtained for other Qm values. 



leading to 83 the number of GRBs in the fiducial HD, i.e. 
20% more data. Moreover, this also makes it possible to re- 
duce the error thanks to the weighted averaging procedure 
employed. Indeed, in the old Schacfer HD, the numbers of 
GRBs with ij,{z) estimated from 1 to 5 correlations are re- 
spectively (3, 12, 14, 27, 13). On the contrary, these numbers 
now read (17, 7, 16, 24, 16, 3) with 3 GRBs having fx com- 
puted from all the six 2D correlations. We do expect that 
this lead to a decrease of the uncertainty on /i(z). Averag- 
ing over the 83 objects, we indeed find fffi/asch = 0.86±0.09 
with {a^/asch) rms — 0.87, i.e. there is on average a signi- 
ficative 14% reduction. 

Summarizing, the detailed analysis performed make us 
confident that the fiducial GRB HD plotted in Fig. 1 and 
tabulated in the Appendix represents a noteworthy step for- 
ward for using GRBs as cosmological tools. 

3.3 Varying the cosmological model 

The fiducial HD has been obtained after calibrating the six 
2D correlations in Table 1 assuming a ACDM concordance 

cosmology. Although this model is in agreement with most of 
the available astrohysical data, there is still not a clear con- 
sensus on the value of its parameters [Qm, h). Moreover, al- 
though not preferred by the data, a varying equation of state 
(EoS) for the dark energy fluid is still a viable (and theoret- 
ically better motivated) option. A large class of dark energy 
models predict a depedence of the EoS on the scale fac- 
tor a which is well fitted by the Ghevallier - Polarski - Linder 
(GPL) ansatz (Ghevallier & Polarski 2001; Linder 2003) : 

w{a) = Wo + Wa{l — a) =wo + ■Waz/{1 + z) . (19) 

The luminosity distance is then computed as : 

rf.(.) = ^(l + .)f ^ (20) 
with E{z) = H{z) /Hq given by : 

E\z) = nM(l + zf + Ox(l + «)3(i+'"o+«.„)g-^ (21) 

where Q.x = 1 — because of the flatness assumption. 

In order to quantify the impact of varying the EoS pa- 
rameters {wo,Wa) and the matter density Q.m, we have re- 
calibrated all the six 2D correlations on a regular grid in 
(wo, Wa) for five different values of Q.m- Namely, we consider 



CPL models with -1.3 < u>o < -0.3 and -1.0 < Wa < 1.0 
for Q,M from 0.20 to 0.40 in steps of 0.05. These ranges are 
larger than what is allowed by the data, but have been cho- 
sen in order to take into account also models very different 
from the fiducial one. For each point in the grid, we repeat 
all the steps described in the previous subsection to get the 
distance modulus to each GRB in the sample. We then col- 
lect these values and evaluate, for each GRB, the root mean 
square of the percentage deviation from the fiducial ^ value. 
The results are plotted in the left panel of Fig. 3. Note that 
in this way we are performing a sort of average of the abso- 
lute percentage deviation so that, for instance, we can say 
that the distance modulus to a given GRB varies of order the 
number in the plot when (wo, Wa) span the range considered 
above. Although averaging over different models is not mo- 
tivated, we have checked that this simple procedure gives us 
the correct order of magnitude of the effect of changing the 
EoS and the matter content on the final Hubble diagram. 

Looking at the left panel of Fig. 3, it is clear that the 
uncertainty on the background cosmological model to use in 
the calibration procedure is quite modest. Indeed, the dis- 
tance modulus may be under or overestimated by a modest 
0.5% with values never larger than 1%. This is is the same 
order of magnitude and often smaller than the measurement 
uncertainty on fi{z) so that one can argue that our fiducial 
GRBs HD may be used even if the underlying cosmological 
model is different from the ACDM one we have assumed. 
Actually, this result is partly due to the use of the distance 
modulus rather than the luminosity distance directly. In- 
deed, as well known, a logarithmic scale always reduce the 
errors, an effect which is at work also in this case as can be 
understood looking at the right panel of Fig. 3. Here, we plot 
the same quantity as the left panel, but referring to the lu- 
minosity distance dhiz). The rms percentage deviation now 
increases to ~ 10% with values as large as 15% so that the 
impact of cosmology is more apparent. However, it is worth 
stressing that it is common practice to use the ijl{z) rather 
than the dL(z) data when constraining cosmological mod- 
els so that we can give off any further consideration of the 
reconstructed luminosity distance. 

It is somewhat surprising that changing the cosmologi- 
cal model used for the calibration has such a small impact 
on the estimated iJ.{z). However, it is easy to trace back the 
origin of this unexpcted result. Figs. 4 and 5 show the depen- 
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Figure 4. The dependence of the a coefficient of the E^-Ep 
correlation normalized to the fiducial value on the Wa parameter 
for different values of wq (namely, too = —1.3, —1.0, —0.7, —0.4 
from top left to bottom right) and Qm (from 0.20 to 0.40 in steps 
of 0.05 from yellow higher to red lower curve). 



dence on the (JIm, wo, ^fa) parameters of the slope coefficient 
a and the intrinsic scatter ai„t (normalized to their fiducial 
values) of the E^y - Ep correlation. Although there are clear 
trends, both a and ai„t change of order few % for the cosmo- 
logical parameters spanning their quite conservative ranges. 
This is also the case for the other 2D correlations we have 
used so that we can deem this as a general result. This re- 
sult then suggests that the calibration is in practice model 
independent, at least within the precision allowed by the 
present measurement uncertainties. As a consequence, the 
final Hubble diagram is almost unaffected by the choice of 
the cosmological model provided this can be well described 
within the framework of the CPL parameterization. 



4 LOCAL REGRESSION CALIBRATION 

Although the above analysis has shown that the choice of 
the underlying cosmological model has only a modest impact 
on the final estimate of the distance modulus, it is worth 
stressing that such a result is far to be definitive. First, we 
have considered only a single (although large) class of dark 
energy models. Moreover, one can argue that future data will 
be affected by smaller measurement uncertainties so that the 
impact of cosmology may turn out to be more important. 
It is therefore of vital importance to look for a method to 
calibrate the 2D correlations through a method that do not 
rely on the assumption of any cosmological model. 

Due to the lack of a set of low redshift GRBs at z < 0.1 
which are cosmology independent, the so called circularity 
problem arise : in order to use the GRBs as cosmological 
tools, one has to calibrate the 2D correlations, but a cosmo- 
logical model has to be assumed for the calibration. In prin- 
ciple, such a problem could be avoided in two ways. First, 
a solid theoretical model must be found in order to phys- 
ically motivate the observed 2D correlations thus setting 
their calibration parameters. Unfortunately, finding such a 
model ramains an ambitious, but still unsuccessful task. Al- 
ternatively, it has been proposed to avoid the need for any 
distance determination by using a sufficiently large sample 
of GRBs within a small redshift bin centred at a whatever 
z (Ghirlanda et al. 2006; Liang & Zhang 2006). However, 



Figure 5. Same as Fig. 4, but for the intrinsic scatter ai^t- Note 
that Qui takes the same values as before, but now increases from 
the yellow lower to the red higher curve in each panel. 



this method might be unrealistic, since the current sam- 
ple of observed GRBs is not large enough. It has recently 
been suggested that such a model independent calibrations 
may be carried out using SNela as distance indicator based 
on the naive observations that a GRBs at redshift z must 
have the same distance modulus of SNela having the same 
redshift. Interpolating therefore the SNela Hubble diagram 
gives therefore the value of n{z) for a subset of the GRBs 
sample with z < lA which can then be used to calibrate 
the 2D correlations (Kodama et al. 2008; Liang et al. 2008; 
Wei & Zhang 2008). Assuming that this calibration is red- 
shift independent, one can then build up the Hubble dia- 
gram at higher redshifts using the calibrated correlations 
for the remaining GRBs in the sample. We present here a 
similar approach still using the SNela as secondary distance 
indicator, but avoiding the need for interpolating the sparse 
dataset through the use of the local regression method we 
briefiy describe in the following. 



4.1 Local regression on SNela 

As a general rule, any interpolation procedure may be seen 
as an attempt of recovering the true underlying curve from 
a sparse dataset. In a sense, interpolating methods try to 
extract the order of a smooth model from the disorder of 
the data originating from that same model. Interpolation 
may therefore be seen as a smoothing procedure so that one 
can naturally resort to the well tested strategies developed 
in this field. A particularly efficient method is known as local 
regression (see Loader 1999 and references therein). 

Originally proposed by Cleveland (1979) and further 
developed by Cleveland and Devlin (1988), the local regres- 
sion technique combines much of the simplicity of linear least 
squares regression with the fiexibility of nonlinear regression. 
The basic idea relies on fitting simple models to localized 
subsets of the data to build up a function that describes 
the deterministic part of the variation in the data, point by 
point. Actually, one is not required to specify a global func- 
tion of any form to fit a model to the data so that there 
is no ambiguity in the choice of the interpolating function. 
Indeed, at each point, a low degree polynomial is fit to a 
subset of the data containing only those points which are 
nearest to the point whose response is being estimated. The 
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Table 2. Calibration parameters (a, 6) and intrinsic scatter aint for the 2D correlations logR = alogQ + b for GRBs with 2 < 1.4 and 
distance modulus estimated from the local regression on SNela. Columns as in Table 1. 

Id M {a,b,a,nt)ML Smed &rms C{S,z) 0.t\l til ('^int)t\l ±1^ 

Ej-Ep 12 (1.53,50.59,0.12) 0.10 0.22 0.22 ^-^^to. It to It 0.19l°;Jj tgj^ 
L-Ep 25 (1.40,52.20,0.40) 0.06 0.45 0.12 l-40j:";2t lo.g? O-^Slalo t'u.'S 



L-Tiag 13 (-0.89,52.07,0.32) 0.11 0.35 0.33 -0.90Io 27 to.64 0-42_o.i3 -0.22 



+0.27 +0.61 p, .r,+0.23 +0.61 
-0.27 -0.64 '-'•*^-0.13 -0.22 

+0.27 +0.57 r, cn+O-lS +0.36 



L-TRT 24 (-0.97, 52.40, 0.44) 0.04 0.47 0.33 
L-V 19 (1.14,52.33,0.54) -0.04 0.56 0.24 ^■'^^to.ll tl.tl O-^^Iq.^s to.2B 



polynomial is fit using weighted least squares with a weight 
function which quickly decreases with the distance from the 
point where the model has to be recovered. 

We apply local regression to estimate the distance mod- 
ulus ij.{z) from the most updated SNela sample. Referred to 
as Union (Kowalski et al. 2008), this compilation includes 
recent large samples from SNLS (Astier et al. 2006) and 
ESSENCE (Wood-Vasey et al. 2007) surveys, older data 
sets and the recently extended data set of distant SNela ob- 
served with HST, all of which have been liomogeriously rean- 
alyzed with the same lightcurve fitter. After selection cuts 
and outliers removal, the final sample contains 307 SNela 
spanning the range 0.015 < z < 1.55. We use this largo 
dataset as input to the local regression estimate of n{z) fol- 
lowing the steps schematically sketched below. 

(i) Set a redshift z where //(z) have to recovered. 

(ii) Order the SNela according to increasing value of 

\z — Zi\ and select the first n = ceJ\fsNeia with a a user 
selected value and A/s Neia the total number of SNela. 

(iii) Define the weight function : 

c {i-\u\'r \u\<i 

W{u) = \ (22) 
I \u\>l 

where u = \z — Zi\/A and A the maximum value of the 
\z — Zi\ over the subset chosen before. 

(iv) Fit a first order polynomial to the data selected at 
step (ii) weighting each SNela with the corresponding value 
of the function W{u) and take the zeroth order term as best 
estimate of ij,{z). 

(v) Estimate the error on ij,{z) as the root mean square 
of the weighted residuals with respect to the best fit zeroth 
order term. 

It is worth stressing that both the choice of the weight func- 
tion and the order of the fitting polynomial are somewhat 
arbitrary. Similarly, the value of a to be used must not be 
too small in order to make up a statistical valuable sam- 
ple, but also not too large to prevent the use of a low order 



polynomial. In order to find what is the better value for 
a and simultaneously check that our choices for the poly- 
nomial degree and weight function do not affect the recon- 
struction, we have performed an extensive set of simula- 
tions. To this aim, we use the CPL EoS and set the model 
parameters (Om, wo, Wa, /i) randomly extracting from the 
conservative ranges 0.15 < Om < 0.45, —1.5 < wo < —0.5, 
-2.0 <Wa< 2.0, 0.60 <h< 0.80. For each redshift value 
in the Union sample, we extract /i(«i) from a Gaussian dis- 
tribution centred on the theoretically predicted value and 
with a standard deviation a = 0.15, in agreement with the 
estimated intrinsic scatter of the SNela absolute magnitude. 
To this value, we then attach an error in such a way that 
the relative uncertainty is the same as the corresponding 
point in the Union sample thus ending up with a mock cat- 
alogue having the same redshift and error distribution of the 
actual one. This mock catalogue is used as input to the rou- 
tine sketched above and finally the reconstructed ij.{z) value 
for each point in the catalog are compared to the input one. 

Defining the percentage deviation Sfi/fi — 1 — 
fj,ir{z)/fj,cPL{z) with fj,ir and hcpl the local regression esti- 
mate and the input CPL values respectively and averaging 
over 250 realizations, we find that the above routine with the 
choice Q = 0.025 gives {5fj,/fi)rms — 0.35% with \ 5fj,/fj,\ < 1% 
independent on the redshift z. Moreover, there is no corre- 
lation at all of Sfi/fj, with any of the cosmological parame- 
ters. We therefore safely conclude that the local regression 
method allows to correctly recover the underlying distance 
modulus at a whatever redshift z from the Union SNela 
sample whichever is the background cosmological model. 

4.2 Model independent calibration 

Having found am efficient way of estimating the distance 
modulus at redshift a in a model independent way, we can 
now calibrate the GRBs 2D correlation using a different ap- 
proach. First, we consider only GRBs with z < 1.4 in order 
to cover the same redshift range spanned by the SNela data. 
As an example, let us consider the L-Ep correlation. For 
each GRBs with measured {z,Pboio), one can write: 

log L = log {4nPtoio) -I- 2 log dz, («) 

= log(47rP6oio) + (2/5)[wr-(«)-25] , 
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Figure 6. The calibrated Hubble diagram with overplotted the distance modulus predicted by the fiducial ACDM model. 



where we have used Eq.(16) and approximated the true 
with its local regression estimate fj,ir{z) without any signi- 
ficative loss of precision. The uncertainty on log L is then 
evaluated as usual by propagation of errors so that we fi- 
nally end with a subset of GRBs with model independent 
estimates of log L that can be used as input to the Bayesian 
fitting procedure described in Sect. 2.2. Carrying on this ex- 
ercise for the first five 2D correlations in Table 1, we get 
the results summarized in Table 2. Note that we cannot 
apply this method to the Lx - Ta correlation since there 
are only 3 GRBs with z < lA after the selection cuts on 
log \Ta/{l + z)] and \ogLx so that the fitting fails. 



Comparing the values in Tables 1 and 2, one should 
get some insight on the dependence of the calibration coef- 
ficients on the redshift. Taken at face values, it is apparent 
that fitting only z < lA GRBs steepens all the 2D correla- 
tions, while the intrinsic scatter is almost unchanged in all 
cases but the - Ep correlation. While this could naively be 
interpreted as an evidence for evolution with redshift of at 
least the slope, such a result loses its statistical significance 
when one looks at the marginalized constraints. Indeed, the 
68 and 95% confidence ranges are quite larger than in the 
previous case so that inferring any constraints on the evolu- 
tion of a is meaningless. Moreover, the present calibration is 
model independent, while the results in Table 1 have been 
obtained assuming a fiducial ACDM model. Although the 
discussion in Sect. 3.3 and Figs. 4 and 5 suggest that this 
should not be a serious problem, it is nevertheless worth 
being cautious before drawing any definitive conclusion. 



4.3 The calibrated Hubble diagram 

We can now use the calibration parameter derived above to 
infer the value of log R (with R = E^ or L) to construct a 
new GRBs Hubble diagram following the same method used 
in Sect. 3.2. The result is what we refer to as the calibrated 
GRBs HD. Plotted in Fig. 6 and tabulated in the Appendix, 
this HD now contains only the 69 Schaefer GRBs since we 
have not used the Lx - Ta correlation being impossible to 
calibrate it with the local regression based method. 

A visual comparison of Figs. 1 and 6 shows that, al- 
though containing a different number of objects (83 vs 69), 
the fiducial and calibrated HDs are in well agreement. This 
is not unexpected since the fiducial ACDM model fits quite 
well the SNela data so that the predicted values for ^{z) 
are in very good agreement with the local regression esti- 
mates based on the same SNela. This qualitative agreement 
is also confirmed by considering the numbers of GRBs that 
are more than 2a away from the theoretical fiducial ACDM 
HD. For the fiducial HD, we find 14 GRBs deviating more 
than 2a, i.e. 17% of the sample, while this number reduces 
to 10, i.e. 14%, for the calibrated HD dataset. Moreover, 
the GRBs are almost the same^ thus arguing that there 



t The likely outliers GRBs for the fiducial dataset are : 990123, 
991208, 991216, 000210, 010222, 020903, 030329, 030528, 050126, 
050406, 060108, 060206, 060614. For the calibrated dataset, the 
list contains : 990123, 990506, 991208, 991216, 000210, 020813, 
050406, 050603, 060108, 060206. There are thus 7 GRBs in com- 
mon within the two lists. Note that 050603 is 3cr away in both 
datasets, while 060614 deviates more than 3a from the fiducial 
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is some problem with the data on this GRBs rather than 
with the procedure used to make up both the fiducial and 

the calibrated HDs. However, since we do not actually know 
whether the adopted fiducial ACDM model is indeed the 
correct one, we prefer to not reject these points in absence 
of a clear alternative motivation. 

To quantitatively compare the fiducial and calibrated 
GRBs HDs, wo have used a linear regression (not taking 
into account the errors on the two variables) to determine 
the slope of the Hfid-^cai relation, with fifid and Hcal the 
distance modulus in the fiducial and calibrated HDs respec- 
tively. We find : 

At/id = (1.02±0.02)Atcai + (-1.28±0.94) 

so that the slope is consistent with 1. It is worth stressing 
that the zcropoint of this relation does not signal a sys- 
tematic offset, but is only a consequence of the slope being 
slightly different from 1. Indeed, defining A^u = /u/id — Mcai 
and the averaging over the 69 GRBs in common, we find 
(A/j) = 0.20 which is consistent with a null value if one 
takes into account the statistical errors on /x/m and Hcai- 

Motivated by this successful comparison, we therefore 
conclude that the local regression based calibration has not 
biased anyway the derivation of the HD so that both the 
fiducial and the calibrated GRBs HDs can hence be used as 
valid tools to constrain cosmological models. 



5 POSSIBLE SYSTEMATICS AND BIASES 

While the previous analysis has shown that the calibrated 
and fiducial HDs arc consistent with each other, such a test 
does not tell us nothing about possible systematic error re- 
lated to the derivation of the underlying GRBs 2D correla- 
tions. Although a detailed analysis of the different system- 
atic errors that can affect these relations is outside our aims 
here, we nevertheless investigate some related issues in or- 
der to qualitatively infer whether they can bias or not the 
derived distance moduli. 



HD only. Wc also stress that only throe GRBs (namely 050603, 
060108, 060206) have ij,{z) estimated also resorting to the Lx - Ta 
correlation so that the use of this new empirical law can not be 
considered as a possible bias. 



5.1 The Lx-Ta correlation 

The use of the Lx - Ta correlation have allowed us to in- 
crease the GRBs sample from 69 to 83 objects^. However, 
since this is the first time such a correlation is used to build 
a GRBs HD, it is worth whehter it is in agreement with the 
previous results. To this end, we have compared the esti- 
mated ^,{z) for the 17 GRBs common to the Schaefer and 
Willingale et al. samples. This is shown in the left panel of 
Fig. 7 where we plot fiDCC vs fisch, i-c. the distance moduli 
as computed from the Lx - Ta and the other 2D correla- 
tions, respectively. Although with a large scatter, the data 
are consistent with a linear relation with slope 1, as it is also 
demonstrated by the right panel where fincc //J-Sch is next to 
1 within 5% for all the points, but two (the only marginally 
deviant cases being GRB 050603 and GRB 060115). Aver- 
aging over the 17 GRBs, we find Hdcc / fJ-Sch = 1.01±0.03 
with {fiDCC /tJ'Sch)rms = 1-01 and no correlation with red- 
shift. As a further test, a linear regression (not weighting 
the points with their errors) gives : 

fJ.DCC = (0.9±0.3)MScft + (5.5±13.5) 

so that the slope is indeed consistent with 1. The small num- 
ber of objects and the large intrinsic scatter makes the er- 
rors on the coefficient extremely large so that the offset de- 
tected is fully devoided of any stastistical meaning. However, 
that such an offset is not present is also suggested by noting 
that (fisch — fJ-DCc) = —0.35 ± 1.3 so that indeed not sta- 
tistically meaningful systematic offset is detected. Repeat- 
ing the same linear fit but excluding the two points with 
fiDCc/fJ-Sch > 1-05 gives: 

fiDCC = (1.2±0.3)^sch + (-8.6±12.2) , 

{fisch - fiDcc) = -0.02 ± 1.03 , 

thus showing that they are not biasing the comparison. 

Wc therefore conclude that the use of the Lx - Ta corre- 
lation does not introduce any bias and hence can be invoked 
to increase the number of GRBs in the fiducial HD. On the 
other hand, should future data disprove the Lx - Ta corre- 
lation, one should reject from the fiducial HD the 14 GRBs 

§ Note that, here, wc only refer to the fiducial GRBs HD since, 
as explained before, the local regression calibration of the Lx - Ta 
correlation is not possible so that this 2D law has not been used 
in the derivation of the calibrated HD. 
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Figure 8. The ratio among the n-i and Hcal values with defined in the text and Hcal the calibrated distance modulus. 



having fi determined by this empirical law only so that one 
gets back the Schaefer Hubble diagram. 

5.2 How many correlations ? 

Although a lot of attention has been devoted recently to 
look for empirical correlations among GRB properties, it is 
still not clear whether the 2D laws we have used here are 
actually all real correlations or a consequence of different 
observational selection effects. Moreover, it is worth not- 
ing that five of them involve a luminosity measurement (the 
bolometric one, L, or the X-,ray one, Lx) and three of them 
rely on a time measurement (being either Tiag, trt or To) 
so that one should also investigate whether they are actu- 
ally independent relations. Investigating this issue in detail 
is outside our aims here, but we can nevertheless estimate 
what should be the impact of giving off one of the relations. 



To this end, we consider the calibrated^ HD and select the 
13 GRBs having p estimated by all the canonical five cor- 
relations. We then denote with /x_i the value of n obtained 

by excluding the correlation i with i = 1, . . . , 5 for the E-, - 
Ep, L-Ep, L-Tiag, L-trt and L-V respectively. We also 
denote by /i-45 the value obtained by excluding the L - trt 
and L-V correlations since they are those with the largest 
intrinsic scatter. Fig. 8 shows that the ratios fi-i / jicai (with 
/icai the value reported in the HD) is very close to 1, with 
only two points deviating more than 4% from unity. More- 
over, there is no apparent differences between the different 
panels and no trend with the redshift or the distance mod- 

^ We prefer to use the calibrated rather than the fiducial HD in 
order to minimize the impact of choosing a cosmological model 
to determine the parameters of each correlation. 
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ulus. Although a more quantitative analysis with a larger 
number of objects is needed to draw a definitive answer, we 
may, however, argue from this prehminary analysis that the 
exclusion of one or two correlations does not alter the re- 
sults. Put in other words, should one of the canonical 2D 
laws employed to get the calibrated HD be dismissed be- 
cause of future evidences, we should not be afraid of any 
bias in the derived distance moduli. 



5.3 The high z GRBs 

Our sample includes 3 GRBs with z > 5, namely 060522 at 
z = 5.11, 050904 at z = 6.29 and 060116 at z = 6.60. As is 
apparent from Fig. 1, there is a clear gap between the two 
z > 6 GRBs and the rest of the sample so that one could 
wonder whether these two objects share the same properties 
of the other GRBs. 

Actually, even if extreme in redshift, both 050904 and 
060522 follow quite well the correlations they are involved 
in. In particular, 050904 does not deviate from five 2D laws, 
while the observed properties of 060522 make this object in 
good agreement with three 2D correlations. As such, we do 
not expect they bias anyway the fiducial HD. Moreover, be- 
ing at such a high redshift, they are not used in the calibra- 
tion of the 2D empirical laws based on the local regression 
analysis. Nevertheless, their distance moduli in the fiducial 
and calibrated HDs are in very good agreement. We there- 
fore conclude that, notwithstanding the doubts on the va- 
lidity of the 2D correlations at very high z, the inclusion of 
these two GRBs does not bias anyway the results. 

5.4 Selection effects 

Since the first works on the existence of correlations among 
GRBs properties allowing to use these objects as standard- 
izeable candles, serious doubts were raised on the actual va- 
lidity of such relations. Although being very luminous, the 
detection of GRBs is a complicated task so that different bi- 
ases may favour the detection of a particular class of GRBs. 
Should this be indeed the case, one could argue that the 
2D empirical laws are not universal or, what is worst, they 
only come out as a consequence of a detection bias (see, e.g., 
Butler et al. 2007 and Shahmoradi & NemirofT 2009). 

One possible example is represented by the shape of the 
GRB spectrum in the high energy region. GRBs more fluent 
at high energies (in their rest frame) at high redshift will be 
more likely to trigger local satellites since their energy in 
redshifted in the satellite bandpass. Approximating the en- 
ergy spectrum as a power -law, i.e. N{E) oc for large 
E, according to this scheme, high z GRBs with more shallow 
profiles will likely be detected preferentially over the steeper 
one. As a consequence, a plot of f3hr vs z should detect no 
correlation even if this is present. Schaefer (2007) reports the 
values of Phr for 28 of the 69 GRBs in the calibrated HD, 
showing that there is indeed no variation with the redshift. 
For instance, roughly averaging over the 28 GRBs with mea- 
sured Phr, we find {l3hr) = —2.42, while averaging over only 
the 17 GRBs with z < 1.4 gives (Phr) = -2.47. Note, how- 
ever, that both distributions are quite asymmetric with most 
of the GRBs having Phr > —2.7 and only 3 with P^r < —3. 
It is also worth noting that, taken at face values, we find a 



clear correlation between log L and Phr with log L increas- 
ing with Phr , although the large errors on Phr prevents from 
deeeming as statistically significant such a trend. 

Quantifying the effect of a selection effect on Phr would 
need the values of this parameter for all the GRBs in the 
sample which is not the case. However, since Phr and log L 
are positevely correlated, one can argue that selecting only 
GRBs with small \Phr\ is the same as selecting those with 
the larger values of log L so that we can try investigating the 
effects of a possible selection on log L as an alternative ap- 
proach. Rather than implementing a procedure based on the 
construction of mock samples (which would also need the de- 
tails of the detection strategies of the different satellites), we 
make an ideal experiment using the GRBs in our sample. As 
a first step, we select only the GRBs with z < 1.4 and divide 
them in three samples referred to as CalAll, Cal51 and Cal52 
containing all the GRBs with logL > A with A = 0,51,52 
respectivelyll . Using these samples, we recalibrate the five 
2D correlations using the local regression method described 
in Sect. 4.2. We now use these calibrated relations to infer 
the GRBs distance moduli whatever the value of log L is 
thus siumlating the case of an unknown selection effect on 
log L being at work. Comparing the distance moduli thus 
obtained, we find : 

HCalSl = {1.00±0m)tICalAll + (-0.1±1.3) , 

f^Cal52 = (0.93±0.06)MCa!AU + (-2.6±2.7) , 
while for the offset we find : 

(MCaiSl - t^CalAll) = 0.29 ± 0.44 , 
{HCal52 - HCalAu) = 0.76 ± 0.91 . 

It is worth stressing that the above linear fit has been done 
without taking into account the errors on the variables thus 
giving all the points the same weight. However, this qualita- 
tive analysis shows that, if a selection on logL should be at 
work, neglecting its effect does not bias the derived distance 
moduli in a statistically meaningful way. On the other hand, 
there is actually an increase in A/j as A increases thus mo- 
tivating further and more quantitative investigation on this 
issue through a careful set of simulations. 



6 GRBS HDS AS A DARK ENERGY PROBE 

The search for standard candles has been always motivated 
by the need of building up reliable Hubble diagrams in or- 
der to probe the evolution of the universe. As Eqs.(16) and 
(20) clearly show, the Hubble diagram allows to estimate 
the luminosity distance as a function of redshift and put 
constraints on the integrated Hubble parameter H{z) and 
hence on the cosmological parameters of a given model. 

The importance of SNela for such a project is impos- 
sible to underrate as it is witnessed by the discovery of the 
cosmic speed up and the overwhelming flood of papers us- 
ing this dataset to test the viability of different cosmological 
models. However, their main limitation is intimately related 

li Actually, tlie CalAII and Cal51 samples are almost identical, 
but we have preferred to not furtherly divide the sample because 
of too low statistics. 
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to their nature preventing us from detecting SNela to red- 
shift higher than z ^ 2. Moreover, even with future satelhte 
missions, the most of the sample wiU be made out of ob- 
jects with z < 1. As such, while they are excellent probes 
over the range (0, 1) following the transition from the dark 
energy dominated present day universe to a decelerating 
expansion, SNela are unable to test the universe during the 
matter dominated era at z > 2. On the contrary, GRBs 
are ideal tools to probe the high redshift universe. That 
GRBs and SNcIa arc complementary may be easily under- 
stood by considering, for instance, the redshift distribution 
of the nowadays sample. Indeed, while SNela mainly cov- 
ers the range (0, 1), GRBs arc observed up to z ^ 6.6 with 
most of them lying at z > 1. As such, one may anticipate 
that they are less sensitive to the details of the dark energy 
EoS since they mainly probe the high z regime deep in the 
matter dominated era. It is this peculiar feature, opposite 
to what happens for SNcIa, that makes them a combined 
SNela -|- GRBs Hubble diagram an ideal tool to constrain 
the universe expansion. 

6.1 Present day data 

In order to show the power of a combined SNela -|- GRBs HD 
as a cosmological probe, we present here a first analysis fit- 
ting this datasct to two popular models, namely the ACDM 
and GPL EoS. For both cases, the scaled Hubble parameter 
E{z) = H{z)/Ho is given by Eq.(21 with the ACDM being 
obtained by setting (wo, Wa) = (—1,0). In order to constrain 
the model parameters, we maximize the following likelihood 
function : 



£(p) oc exp (-X /2) , 



(23) 



with 
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In Eq.(24), //ots and ^ith are the observed and theoretically 
predicted values of the distance modulus given the set of 
model parameters p, while the sum is over the MsNeia = 307 
SNela in the Union dataset and the Mgrb = 69 GRBs in the 
calibrated HD. Note that, for GRBs, we have added to the 
denominator the constant error parameter Oint to take into 
account the intrinsic scatter** around the Hubble diagram 
due to the scatter in the 2D correlations. The last two terms 
in Eq.(24) axe Gaussian priors on h and ujm = 0,Mh? and 
are included in order to help breaking the degeneracy with 
the other model parameters. To this aim, wo resort to the 
result of the HST Key project (Freedman at al. 2001) to set 
(hnsT, cthst) = (0.72, 0.08) and to the WMAP5 constraints 

** Note that a similar term is also present for SNela, but it is 
estimated to be (Tint = 0.15 and yet included into the error ai 
provided in the dataset. 



(Dunkley et al. 2008) giving (w^^^^, ctm) = (0.138, 0.004). 
Note that we have not added further data from other probes 
since we are here more interested in testing the ability of the 
combined HD to probe models rather than setting strong 
constraints on the model themselves (which will be the sub- 
ject of a forthcoming publication). 

In order to efficiently explore the parameter space, we 
run a Markov Chain Monte Carlo code generating chains 
with more than 100000 points which are enough to reach 
convergence. We then cut out the first 30% of the chains to 
avoid the burn in phase and thineed the chains taking one 
point each ten points in order to reduce fake correlations. 
For the ACDM model, we find for the best fit parameters : 

inM,h,(Tint) = (0.283,0.701,0.44) 

giving 

XSNela/d.O.f. = 1.03 , XGRs/d.O.f. = 1.11 , COM = 0.139 

with d.o.f. = JVdata - Mp and ATp = 3 for the ACDM case. 
The reduced close to 1 for both the SNela and GRB 
datasets and the very good agreement with the WMAP5 
ujm value clearly shows that the concordance ACDM model 
still stands out as an excellent description of the universe 
expansion. After marginalizing over the other quantities, the 
constraints on the model parameters are found to be : 



n 007+0. 013 +0.027 
U.ZO(_o _o.029 



h = 0.700 



+0.006 +0.012 
-0.006 -0.012 



Cint 



f. ir+0.21 +0.40 
U.10_Q _0.13 



where, following the Bayesian prescription, we give the me- 
dian as central value and the reported errors refer to the 
68 and 95% CL. It is worth noting that the constraints on 
aint shows that this quantity is weakly constrained with the 
median values being very different from the best fit one. 
This is, however, not surprising. On the one hand, it is well 
known that, because of degeneracies and projection effects, 
the maximum of jC{p) is not given by the median values of 
the parameters. This could be the case only if the covariance 
matrix of the parameters is diagonal which is quite unusual. 
On the other hand, we have not included in the likelihood 
definition any prior on aint which could be inferred consider- 
ing the intrinsic scatter of the 2D correlations used to build 
the GRB HD. Moreover, the asymmetry in the aint CL is a 
result of the obvious constraint that this quantity must be 
positive because of its same defintion. 

As a further example, we now consider the CPL case 
leaving the parameters (uio,Wa) free in the fit thus giving 
J\fp = 5. For the best fit, we find : 

{Qm, wo, Wa,h.ai„t) = (0.273, -1.28, 2.02, 0.707, 0.45) 

giving 

XSNela/d.O.f. = 1.03 , XGRs/d.o.f. = 1.06 , Wm = 0.136 

so that there is still a perfect agreement with the data. Note 
that, formally, the CPL model fits the GRB data better than 
the ACDM one, but the increase by two of the number of 
parameters is actually not compensated by a corresponding 
reduction of the reduced Xgrb so that a pure cosmological 
constant is still preferred. One could be surprised that the 
best fit parameters seem to favour a model which is very 
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different from the ACDM one with the present day value 
of the EoS deep in the phantom regime wo < — 1 and a 
strongly varying EoS. However, this is just a consequence 
of the parameter degeneracies quoted above^^. Indeed, the 
Bayesian constraints turn out to be : 

O — n 900+0 021 +0.036 
i'M — '-'.^i50_Q 022 -0.039 1 



measurement. For SNela, following Kim et al. (2004), we 
adopt'^'l" : 



Wo 



, nfi+019 +0.34 
-l.UD_0 2l -0.42 5 



Wa — U.O_i 1 _i 9 , 

^ _ n fioq+o oie +0.029 

a — U.U3tl_Q _o.022 ; 



n ns+015 +0.42 

■ '-'•'-'"-O.OS -0.07 I 



showing that the ACDM model with {wo, Wa) = (—1,0) is 
well within the 68% CL. It is worth noting that the con- 
straints ou both iL'o and Wa are quite weak, as those on the 
intrinsic scatter aint- This same result is obtained if one 
only uses the SNela HD so that one could wonder whether 
adding the GRB data have improved the situation. Actu- 
ally, this is not becuasc, as we have yet quoted above, the 
GRBs HD only probes the matter dominated era when the 
dark energy term is almost disappeared so that the pre- 
cise values of {wo,Wa) are meaningless. As we will see later, 
GRBs improve, on the other hand, the constraints on Qm 
even if this is not apparent for the fits we are considering 
here because of the addition of the priors on h and ujm- 
As a final remark, we note that, even with these weaJi con- 
straints, we can exclude with high confidence historically 
discussed cosmological models with, e.g., (f2M,f^A) = (1,0) 
or {Qmj^a) = (0.3,0), and have a further strong evidence 
for an accelerating universe. 



6.2 Future SNela -h GRB data 

The analysis of the present day data have not shown the 
full complemetarity between SNela and GRBs, essentially 
because we have added the priors on the two parameters h 
and ujm- To better investigate this issue, we compute the 
Fisher information matrix (see, e.g., Tegmark et al. 1997 
and references therein) defined as : 



dpidpj 



(25) 



with L = — 21n>C(p), and we redefine the likelihood as in 
Eq.(23), but excluding the two priors on (/i, wm)- Note that 
(. . .) denotes the expectation value, but, actually, this is 
computed by evaluating the Fisher matrix elements for fidu- 
cial values of p. The inverse of the Fisher matrix gives the 
covariance matrix C so that cr, = \/Cu is an estimate of the 
error on the parameter i. We remember that, because of the 
Cramer - Rao theorem, this is an upper limit on the error 
on the parameter pi that can be obtained by an experiment 
with the given characteristics. 

A key ingredient in the computation of F is represented 
by the details of the SNela and GRB surveys giving the 
redshift distribution of the samples and the errors on each 



Note also that we have not set here the usual prior wo+Wa = 
which is typically introduced to assure that the dark energy EoS 
fades away in the early universe. 



a{z) = \ cr'iya + 



with z,nax the maximum redshift of the survey, asys an irre- 
ducible scatter in the SNela distance modulus and am to be 
assigned depending on the photometric accuracy. For GRBs, 
we assume that the redshift distribution is provided is the 
same as the star formation history in Porciani & Madau 
(2001) although it is worth noting that this just a first order 
approximation. Indeed, one should convolve this theoreti- 
cal distribution with the detection efficiency of a satellite, 
like Swift or Fermi, which will also claim for a model of 
GRB explosion to compute the expected fluence. Since all 
these quantities arc up to know largely uncertain, we prefer 
to not detail them warning the reader that our results are 
somewhat optimistic. Concerning the error distribution, we 
assume that the main source of error is the intrinsic scat- 
ter (as it is indeed the case with the present day data) so 
that we give to each GRBs distance modulus an uncertainty 
randomly generated from a Gaussian centred on acRB and 
with standard deviation (Jstd- To further concentrate our 
attention on the constraints on the cosmological parame- 
ters, we moreover assume that the intrinsic scatter aint in 
Eq.(24) has been estimated somewhat and included in the 
error on fj, so that this is no more a quantity to be con- 
strained. Finally, in order to run the Fisher matrix calcu- 
lation, we have to choose a fiducial model which we do by 
setting {flM,wo,Wa,h) = (0.277,-1,0,0.72) in agreement 
with the WMAP5 results. 

As a first test, we consider the present day data thus 
using the observed redshift distribution as the Union one 
for SNela and that of the calibrated HD for GRBs. Wc fur- 
ther set {MsNeia,am) = (307,0.33) and {aGRB,asTD) = 
(0.58, 0.22). With these values, we get : 

aiClM) = 0.58 , (t(wo) = 0.17 , ct(w„) = 1.2 , 

using SNela data only, while adding GRBs gives : 

a{QM) = 0.08 , (t{wo) = 0.08 , ff(w„) = 0.4 . 

Such a test immediately shows the complementarity of 
SNela and GRBs with these latter data halving the error 
on Qm thanks to their ability of probing the matter domi- 
nated era. Note also the reduction of the errors on {wq, Wa) 
which seems to contradict our claim that GRBs are not sen- 
sitive to the dark energy EoS. Actually, the reduction is due 
mainly to GRBs increasing the statistics rather than probing 
the EoS as is better understood considering future data. 

To this end, we now simulate a SNAP - like SNela survey 
setting {AfsNeia, am) — (2000,0.02) and using the Porciani 
& Madau (2001) distribution to generate 200 GRBs leaving 
the error parameters {aoRBjCrstd) unchanged. Using only 
SNAP SNela, we get : 

a{QM) = 0.077 , cr(wo) = 0.039 , cT{wa) = 0.22 , 



+i Note tliiit. in Kim et al. (2004), the authors assume the data 
are separated in redshift bins so that the error becomes = 
'^lys/^bin + ■N'bin{z/zmax)'^a'^ with jVi,j„ the number of SNela 
in a bin. However, we prefer to not bin the data so that = 1. 
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while adding GRBs gives : 

a{QM) = 0.019 , a(w„) = 0.036 , a{wa) = 0.20 . 

The role of the complementarity of SNela and GRBs is now 
better explained. Indeed, while a{QM) is reduced by a factor 
~ 4, both a{wo) and a{wa) are almost unaltered by the 
addition of GRBs in agreement with the claim that their 
Hubble diagram is unable to put constraints on {wQ,Wa)- 

As a final remark, it is worth noting that the detection 
of 200 GRBs is likely well within the capabilities of both 
the Swift and Fermi satellites so that one can argue that, by 
the time the SNAP survey will be completed, one indeed has 
the opportunity to build up a combined SNela + GRB HD as 
the one assumed in our Fisher matrix forecast. However, it is 
worth stressing that the precision attainable with GRBs may 
critically depend on their assumed redshift distribution and 
on the intrinsic scatter of the 2D correlations used to build 
up the calibrated HD. Both these key ingredients have been 
modelled very roughly in our analysis so that a more careful 
investigation should be needed. Although such a work can 
downgrade our optimistic forecast, we argue that similar 
results can still be obtained by increasing the number of 
GRBs still remaining within the realm of the Fermi satellite 
detection capabilities. 



7 CONCLUSIONS 

As soon as the need to probe the Hubble diagram to a 
redshift higher than the one attainable with SNela be- 
came clear, most attention has been devoted to search for 
a method to make GRBs standardizeable candles. To this 
aim, different 2D correlations between a luminosity or en- 
ergy distance dependent quantity and observationally acces- 
sible GRBs properties have been proposed in order to build 
up a GRBs Hubble diagram to be used for constraining cos- 
mological parameters. In an attempt to make a step forward 
along this road, we have here reconsidered different aspects 
of the problem trying with the final aim of constructing a 
more reliable GRBs Hubble diagram. 

First, we have added the recently found Lx - Ta correla- 
tion to the other five 2D empirical laws used up to now. Al- 
though being different from the other ones since it is based 
on afterglow rather than peak quantities, we have shown 
that the estimated distance modulus for each of the GRBs 
common to the Schaefer (2007) and WilUngale et al. (2007) 
samples is in agreement with the one obtained using the 
standard set of 2D correlations. As such, one can be confi- 
dent that no systematic bias is introduced by resorting to 
this new empirical law. On the contrary, its use makes it 
possible to both reduce the errors on fi{z) by a significa- 
tive ~ 14% and increase the sample from 69 to 83 GRBs. 
With the Swift satellite still operating and collecting data 
on the X - ray afterglow, the prospects to both improve the 
estimate of the calibration parameters and adding still more 
GRBs to the sample are quite good. Needless to say, this is 
also possible for the peak based 2D correlations. It is worth 
noting that, in order to reduce the uncertainty on fi{z), fol- 
lowing Schaefer (2007), we have made a weighted average 
of the estimates coming from the use of the six 2D corre- 
lations. While this makes it possible to, in a sense, average 
out any systematics intrinsic to each particular correlations. 



such a procedure is not statistically well motivated. Actu- 
ally, comparing the estimated n{z) with each other and with 
the weighted average makes us confident that this method 
have not biased anyway the final results. However, increas- 
ing the sample will allow us to bettor cross check the results 
thus leading further confidence in this approach. 

As a preliminary step in the construction of the GRBs 
HD, wc have updated the calibration of the six 2D correla- 
tions we have used. To this aim, we have first to assume a 
background cosmological model to compute the distance de- 
pendent quantities. The recent WMAP5 data analysis has 
motivated our choice of a fiducial flat ACDM model with 
(Q.M,h) = (0.291,0.697), but we have also explored the ef- 
fect of varying both {^lM,h) and the {wo,Wa) parameters 
of the phenomenological CPL dark energy EoS. These data 
are then used as input to a Bayesian motivated fitting pro- 
cedure which makes it possible to estimate both the slope a 
and the intrinsic scatter aint of the 2D correlations correctly 
taking into account the errors on both involved variables. To 
our knowledge, this is the first time that the 2D GRBs em- 
pirical laws have been calibrated using the same statistical 
tool for all of them. Moreover, our method avoids to un- 
derestimate the intrinsic scatter which plays an important 
role not only in the determination of the distance modulus 
uncertainty, but also in testing theoretical models of GRBs 
explosion and afterglow. As a somewhat surprising result, we 
find that the calibration parameters (a, b, aint) very weakly 
depend on the assumed cosmological model, i.e. on the value 
of (JIm, ft, Mo, Wa). As a consequence, the fiducial GRBs HD 
(i.e., the one relying on the calibrations made in the fidu- 
cial ACDM model) is essentially unaffected by varying the 
cosmology. Although this nice finding argues in favour of 
the use of the fiducial HD as a tool to constrain different 
cosmological models, we caution the reader that we have 
only explored a single class of dark energy theories. Indeed, 
the CPL parametrization makes it possible to mimic a quite 
large set of models, but may fail to mimic scenarios with a 
significative quantity of dark energy at high redshift. It is 
likely that the impact on {a, b, aint) is laxger in this class 
of theories, even if it is worth stressing that they are dis- 
favoured by the present dataset and the constraints from 
primordial nucleosynthesis. 

In order to fully avoid the problem of assuming any 
background cosmology in the calibration procedure, we have 
illustrated a novel method relying on the estimate of the 
distance modulus of a given GRBs at redshift z through a 
local regression analysis of the Union SNela sample. As we 
have convincingly shown, local regression allows to recover 
the right value of fi{z) whatever are z and the parameters 
(JIm, h, wo,Wa) of the underlying CPL cosmology. As a con- 
sequence, we can calibrate the 2D correlations using GRBs 
with z < 1.4 (i.e., overlapping the redshift range probed by 
SNela) without the need of assuming any fiducial cosmo- 
logical model. This procedure finally leads us to what we 
term the calibrated GRBs HD containing 69 objects with 
model independent estimates of the distance modulus. Re- 
assuringly, for the GRBs in common, the calibrated and fidu- 
cial HDs are in very good agreement thus suggesting that 
both of them arc not affected by any systematic bias induced 
by the different cafibration procedures. 

One cam wonder which is the GRBs HD better suited as 
a cosmological tool. This is somewhat a matter of personal 
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taste. Indeed, the two HDs agree in all of their properties 
also providing the same percentage of likely outliers with 
respect to a fiducial ACDM model. On one hand, the fidu- 
cial HD contains a larger number of objects with smaller 
uncertainties, but, although this is unlikely, one can still not 
exclude a bias originating from the model dependent cali- 
bration. On the other hand, the calibrated HD is free of any 
problem due to the unknown background cosmology, but its 
constraining power is likely to be weaker because of the lower 
number of objects and the greater uncertainties, both effects 
due to not having used the Lx - Ta correlation. Should fu- 
ture data make it possible to calibrate also the Lx - Ta law 
with the local regression based method, the calibrated HD 
should become as precise as the fiducial HD. As a concluding 
remark, we therefore argue in favour of this latter approach 
as the best way to use GRBs as cosmological tools. 
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APPENDIX A: HUBBLE DIAGRAMS DATA 

As explained in the main text, the fiducial and calibrated 

HDs contain a different number of GRBs because it is not 
possible to use the local regression calibration for the Lx - Ta 
relation. Therefore, we report in Table Al the redshift and 
distance moduli of the 69 GRBs common to both samples, 
while Table A2 gives the same data for the 14 GRBs present 
only in the fiducial HD. Further material (such as the values 
of II from the single 2D correlations) is available on request. 



APPENDIX B: COMMENTS ON OUTLIERS 

As we have quoted in Sect. 4.3, there are 13 (10) GRBs 
deviating more than 2a from the fiducial (calibrated) HD 
which are therefore likely outliers. Actually, dealing with 
likely outliers is a complicated task since one has to find 
a serious motivation before rejecting a given object from a 
sample. To this end, one should look both at its eventually 
peculiar features and the instrumental and observing setup 
which has been used to gather the corresponding data. Since 
we have no access to these informations for all GRBs, we 
have preferred to not reject them. 

Nonetheless, as an example, we discuss here the case 
of GRB 060614 which is the only GRB deviating more than 
3fT from both the fiducial HD and the Lx - Ta correlation. 
Indeed, this object is somewhat very peculiar being the first 
long GRB observed at low redshift (z = 0.125) and clearly 
not associated to a bright Ib/c SN (Delia Valle et al. 2006). 
Actually, its short duration (Tgo = 100 s) makes it hardly 
classifiable as a long GRB, while its morphology is quite 
similar to the one of GRB 0507024 which is indeed classified 
as a short GRB (Piro 2005; Zhang ct al. 2007). It is worth 
noting that GRB 060614 has been interpreted both in the 
framework of the canonical fireshell model (Ruffini et al. 
2001; Caito at al. 2008) and by (Norris & Bonnell 2006) as 



Table Al. Fiducial and calibrated Hubble diagrams data for 
the 69 GRBs common to both samples. Order of the columns as 
follows: 1. GRBs id; 2. redshift; 3., 4. distance moduli (with Icr 
error) from the fiducial and calibrated HDs. 
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Table Al - continued 
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Table A2. Hubble diagram data for the 14 GRBs present in the 
fiducial one only. Order of the columns as in Table Al. 
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a new class of intermediate sources. All these peculiarities 

indeed argue in favour of the rejection of this GRB from the 
fiducial HD as a physical outlier. 

Another interesting case is represented by GRB 060505 
which is not associated by any SN thus suggesting it is rep- 
resentative of a class of objects with different progenitors. 
However, among the likely outliers, we also find GRB 030329 
wich is one of spectroscopically confirmed GRBs associated 
with a SN. However, its X-ray afterglow was observed with 
only two pointings of the Rossi - XTE instrument obtained 
respectivly 5 hours and 1.24 days after the burst (Mar- 
shall & Swank 2003; Marshall et al. 2003). Actually, an in- 
complete coverage of the X - ray lightcurve typically takes 
place for many GRBs observed with the old missions such 
as Beppo SAX, Rossi, Hete2. One can argue that the pa- 
rameters obtained in these cases are somewhat uncorrectly 
estimated thus altering the position of the corresponding 
GRB in the fiducial HD. This could be the case for 10 out 
of 13 likely ouliers GRBs (namely: 990123, 991208, 991216, 
990506, 000210, 010222, 020813, 020903, 030329, 030528) 
thus suggesting that they could be rejected because of ob- 
servationally (rather than physically) motivated problems. 

As it should be clear from these comments, we have up 



to now only some phenomenological and qualitative hints 
to address the problem of outliers. We therefore deserve a 
detailed analysis also taking into account more quantitative 
features to a forthcoming publication. 
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